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ABSTRACT 

The observational frontiers for the detection of high-redshift galaxies have recently 
been pushed to unimaginable distances with the record-holding Lyman Alpha Emitter 
(LAE) UDFy-38f 35539 discovered at redshift z = 8.6. However, the physical nature 
and the implications of this discovery have yet to be assessed. By selecting galaxies 
with observed luminosities similar to UDFy-38135539 in state-of-the-art cosmological 
simulations tuned to reproduce the large scale properties of LAEs, we bracket the 
physical nature of UDFy-38135539: it has a star formation rate ~ 2.7 — 3. 7M© yr _1 , it 
contains 10 8 - 3 - 8 - 7 M Q of stars 50-80 Myr old, with stellar metallicity ~ 0.03 - 0.12Z Q . 
For any of the simulated galaxies to be visible as a LAE in the observed range, the 
intergalactic neutral hydrogen fraction at z = 8.6 must be \hi ^ 0.2 and extra 
ionizing radiation from sources clustered around UDFy-38135539 is necessary. Finally, 
we predict that there is a 70% (15%) probability of detecting at least 1 such source 
from JWST (HST/WFC3) observations in a physical radius ~ 0.4 Mpc around UDFy- 
38135539. 

Key words: cosmology: theory, galaxies: individual, galaxies: high redshift; galaxies: 
intergalactic medium 



1 INTRODUCTION 

The earliest galaxies which ushered in the era of the cos- 
mic dawn changed the state of the intergalactic medium 
(IGM) out of which they formed in numerous ways: they 
polluted it with heavy elements, heated and (re)ionized it, 
thereby affecting the evolution of all subsequent generations 
of galaxies. We are now in the golden era for the search 
for such galaxies, made possible by the use of ingenious se- 
lection techniques. The first of these, the standard dropout 
technique (e.g. Steidel et al. 1996; Giavalisco et al. 2004) 
relies on broad band filters to detect the Lyman break at 
912 A in the galaxy rest frame. Although galaxy candidates 
have been detected upto z ~ 10 (Bouwens et al. 2009) us- 
ing this method, it has the drawback that the exact source 
redshift cannot be determined with complete confidence. In 
this respect, searches for the Lyman Alpha (Lya) line at 
1216 A (in the galaxy rest frame) using narrow band filters 
have been far more successful in confirming the detection 
of high-redshift galaxies; specific spectral signatures includ- 
ing the strength, width and continuum break bluewards of 
the Lya line make the detection of LAEs (galaxies show- 
ing a strong Lya line) rather unambiguous. Narrow band 
searches have enabled the confirmation of hundreds of LAEs 
in a wide high-redshift range, at z ~ 5.7 (Malhotra et al. 
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2005; Shimasaku et al. 2006), z « 6.6 (e.g. Kashikawa set al. 
2006) and z ~ 7 (e.g. Vanzella et al. 2010). Using the same 
technique, recently, a LAE, designated UDFy-38135539, has 
been confirmed at z = 8.6 (Lehnert et al. 2010), making it 
the farthest astrophysical object known so far; it has over- 
taken the redshift record of z = 8.2 set by the Gamma Ray 
Burst GRB090423 (Salvaterra et al 2009; Tanvir et al. 2009). 



This object has already been observed by a number of 
groups: Lehnert et al. (2010) have observed both the Lya 
and ultraviolet (UV) luminosity for UDFy-38135539, and 
Finkelstein et al. (2010) have obtained broad band UV and 
Spitzer data points for this galaxy, designated ID 125 in 
their work. However, the physical nature of the galaxy and 
the cosmological implications of its discovery have yet to 
be assessed. In this work, we use state-of-the-art cosmo- 
logical simulations coupled to a previously developed LAE 
model (Dayal, Ferrara & Gallerani 2008; Dayal et al. 2009; 
Dayal, Ferrara & Saro 2010), tuned to reproduce a number 
of observables of LAEs, to bracket the physical properties 
of UDFy-38135539. We use the observed Lya luminosity to 
get a hint on the ionization state of the IGM at z ~ 8.6, 
calculate how much this galaxy could have contributed to 
reionization and make predictions for the number of Lyman 
Break Galaxies (LBGs) that could be detected in its vicinity. 
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2 THEORETICAL MODEL 

We start from the analysis of a z = 8.6 snapshot of a set of 
cosmological simulations carried out using the TreePM-SPH 
code GADGET-2 (Springel 2005) with the implementation of 
chemodynamics as described in Tornatore et al. (2007). The 
periodic simulation box has a comoving size of 75ft -1 Mpc 
and initially contains 512 3 particles each of dark matter 
(DM) and gas. The masses of the DM and gas particles are 
moM ^ 1.7 x 10 8 ft -1 M o and m gas ~ 4.1 x 10 7 fr -1 M Q , 
respectively. For each "bona-fide" galaxy (which has at 
least 20 bound particles; see Saro et al. 2006 for details) 
in the simulated volume at z = 8.6, we compute the DM 
halo/stellar/gas mass, star formation rate (SFR), mass- 
weighted age, gas/stellar metallicity, and gas temperature. 
The adopted cosmological model for this work corresponds 
to the ACDM Universe with Q m = 0.26, fl A = 0.74, Q b = 
0.0413, n s = 0.95, H = (100ft) = 73 km s -1 Mpc -1 and 
as = 0.8, thus consistent with the 5-year analysis of the 
WMAP data (Komatsu et al. 2009). Complete details of 
these simulation runs can be found in Dayal et al. (2009). 

2.1 Intrinsic luminosities and dust 

Star forming galaxies produce their intrinsic Lya line and 
UV continuum luminosities (L« nt and L'J 1 *, at restframc 
wavelengths X a = 1216 A and A c = 1700 A, respectively) 
both via stellar (and nebular) emission, and cooling radi- 
ation from collisionally excited neutral hydrogen (Hi) in 
their interstellar medium (ISM). While the spectral energy 
distributions (SEDs) are obtained using the population syn- 
thesis code STARBURST99 (Leitherer et al. 1999) to calcu- 
late the stellar and nebular emission, the latter depends on 
the temperature distribution in the ISM gas. The interested 
reader is referred to previous work (Dayal, Ferrara & Saro 
2010) for a comprehensive discussion. 

As Lya and continuum photons are efficiently absorbed 
by dust grains, we calculate the dust content of each galaxy 
by considering SNII to be the main dust producers; this 
is justified by the fact that the typical evolutionary time- 
scale of evolved stars (^ 1 Gyr) becomes longer than the 
age of the Universe above z £ 5.7 (Todini & Ferrara 2001). 
We further pose that: (i) 0.5 Mq of dust is produced per 
SNII (Todini & Ferrara 2001; Nozawa et al. 2007), (ii) SNII 
destroy dust in forward shocks with an efficiency of about 
40% (McKee 1998; Seab & Shull 1983), and (iii) a homo- 
geneous mixture of gas and dust is assimilated into further 
star formation. To calculate the dust optical depth, r, to 
UV continuum photons, based on the results obtained at 
z ~ 6.6, we assume that dust is made up of carbonaceous 
grains and spatially distributed as the stars (Dayal, Ferrara 
& Saro 2010). The fraction of continuum photons escaping 
the galaxy undamped by dust is then f c — (1 — e -T )r -1 . The 
analogous quantity for Lya photons is taken to be f a = 1 so 
as to obtain the maximum possible Lya luminosity emerging 
from the galaxy itself. 

2.2 IGM transmission and source clustering 

After escaping out of the galactic environment, Lya pho- 
tons are further attenuated as they travel through the IGM. 
Depending on the intergalactic hydrogen ionization state, 



only a fraction < T a < 1 of the Lya photons emerg- 
ing out of a galaxy undamped by dust reaches the ob- 
server (continuum photons are instead unaffected by the 
IGM); T a hence depends upon the IGM ionization state. 
Since this value is largely unconstrained, we explore 6 dif- 
ferent values of the average IGM H i fraction at z = 8.6: 
Xhi = 0.2,0.3,0.4,0.5,0.6,0.7. For each value of xhi, wc 
compute the nominal radius, Ri, of the spherical, ionized 
HII region around each simulated galaxy. In reality, though, 
due to source clustering (i.e. when the separation between 
any two galaxies becomes smaller than either of their ion- 
ized region radii), multiple galaxies can contribute ionizing 
photons to the same ionized region, which will then be char- 
acterized by an effective radius, R e / > Ri, calculated as 
follows. Suppose galaxies 'A' and 'B', and 'B' and 'C' have 
overlapping ionized regions. Then, the effective ionized vol- 
ume 'A' is embedded in is the sum of the ionized volumes 
of 'A', 'B' and 'C; the same holds true for both 'B' and 'C'. 
Within this ionized volume, the total photoionization rate 
at distance r from 'A', is 
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where the terms on the right hand side represent the pho- 
toionization rate from (i) the direct radiation from 'A', (ii) 
the galaxies clustered around 'A' and (iii) the ultravio- 
let background (UVB) from distant galaxies, respectively. 
The UVB photoionization rate is related to \hi by Yb — 
(1 — XHi) 2 X~H 1 i nH ° LB w h ere n H is the mean IGM hydrogen 
number density and as is the case-B hydrogen recombi- 
nation coefficient. Further, L e ™ A — L 1 ™ 1 f e3C , is the specific 
ionizing luminosity emerging from 'A' and / eS c = 0.02 is 
the escape fraction of H i ionizing photons (Gnedin et al. 
2008), L e J™ is the ionizing luminosity emerging from the i th 
galaxy of the total 'N' galaxies with which 'A' shares an ion- 
ized region, vl is the frequency corresponding to the Lyman 
limit wavelength (912 A), gl is the hydrogen photoioniza- 
tion cross-section and Ti A is the distance between galaxies i 
and 'A'. This procedure is carried out for each galaxy in the 
simulated volume. We then assume photoionization equilib- 
rium to compute \hi within the effective ionized region of 
each galaxy. At the edge of this region, we force \hi to attain 
the assigned global value. We use the complete Voigt profile 
to calculate the optical depth of Lya photons along the line 
of sight to compute T a . The observed Lya and UV contin- 
uum luminosity are then simply L a = L™ f a T a = L™ T a 
(the latter equality descends from our maximal assumption 
f a = 1), and L c = L l c nt f c , respectively. 

As a final step, we select galaxies that would be ob- 
servable as LAEs according to the canonical criterion, L a 
10 42 ' 2 ergs -1 and observed equivalent width EW ^ 20 A. 
Among these, we further isolate those that fall within the 
La and continuum magnitude Mjjv range observed by Lehn- 
ert et al. (2010) and Finkelstein et al. (2010): using SIN- 
FONI, the former find L a = (2.7 - 8.3) x 10 42 ergs -1 , 
M uv = -19.6 to -18.6, while Finkelstein et al. (2010) find a 
tighter bound of Muv = —19.2 to —19.0. To summarize, all 
simulated galaxies that have L a — (2.7 — 8.3) x 10 42 ergs -1 
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Figure 1. The UV magnitude plotted as a function of the ob- 
served Lye? luminosity for LAEs at z = 8.6. The yellow shaded 
area delimitates the region of the parameters deduced from the 
observation of UDFy-38135539 by Lehnert et al. (2010), the green 
shaded area shows the UV magnitude range measured for the 
same galaxy by Finkelstein et al. (2010) and points represent 
galaxies identified as LAEs in our simulation; the 7 LAEudf for 
Xhi = 0-2 are marked by triangles. Arrows show the effects of 
different physical processes on the observed luminosity (see text 
for details). 
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-19.2 to —19.0 and designated as LAEudf in 



3 RESULTS 

We now discuss the main results obtained from the above 
mentioned computations. We find that, without requiring 
any tuning of the model parameters, no LAEudf are found 
in the simulated volume for xhi > 0.2. For xhi = 0.2, we 
find a total of 215 LAEs in the simulated volume, 7 of which 
are identified as LAEudf as shown in Fig. 1. From now 
on, LAEudf refer to the LAEs in the combined observed 
Lya and UV luminosity range of Lehnert et al. (2010) and 
Finkelstein et al. (2010) for xhi = 0.2. 

We digress here to discuss the two main ingredients 
whose combined effects can change the slope of the L a — 
Muv relation shown in Fig. 1. The first of these concerns the 
dust attenuation: for any given LAE, an increase (decrease) 
in f c leads to the galaxy becoming brighter (fainter) in the 
UV, shifting the position of the galaxy vertically upward 
(downward) on the plot; a decrease in the value of f a (recall 
that we have used f a — 1 in our model) due to dust atten- 
uation of Lya photons leads to a corresponding decrease in 
L a , moving the points horizontally leftward on the plot. The 
second effect is the change in T a due to peculiar velocities: 
inflows (outflows) of neutral gas into (from) the galaxy im- 
part an extra blueshift (redshift) to the Lya photons, lead- 
ing to a decrease (increase) in f a (Santos 2004; Verhamme 
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Figure 2. Comparison of the theoretical SEDs (lines) for the 7 
LAEudf compared to the data (points) collected by Finkelstein 
et al. (2010); the points at 3.6 and 4.5^m have been obtained by 
Spitzer observations and downward pointing arrows show the la- 
upper limits. 



et al. 2006; Dijkstra & Wyithe 2010; Dayal, Maselli & Fer- 
rara 2011) moving the relation horizontally towards the left 
(right). 

Further constraints on the LAEltdf come from the 
SEDs which we have obtained from STARBURST99 includ- 
ing dust attenuation (f c ranges between 0.24-0.3 for the 
LAEudf) using the supernova extinction curve (Bianchi & 
Schneider 2007). As shown in Fig. 2, the model SEDs are 
in excellent agreement with the observed data points, in- 
cluding the two from Spitzer at 3.6 and 4.5/xm obtained by 
Finkelstein et al. (2010). The SED is a crucial constraint on 
the physical properties of galaxies that could be identified 
as LAEudf' although for values of xhi < 0.2 (for the given 
dust model) a larger number of galaxies would fulfill the 
LAEudf selection criterion, their physical properties can- 
not vary too much without the SEDs becoming inconsistent 
with the observed one. 

As the 7 LAEudf galaxies closely resemble the ob- 
served properties of UDFy-38135539, it is reasonable to 
conclude that their physical properties should also match 
those of that object, which are now discussed. LAEudf 
show SFRs in the range 2.7 — 3.7M0yr _1 , as seen from 
Panel (a) of Fig. 3. These values are quite comparable to 
the ones inferred by Lehnert et al. (2010) from the Lya lu- 
minosity of UDFy-38135539 (0.3 - 2.1 M yr" 1 ). LAEudf 
have stellar masses ~ lO 8 3-8 ' 7 M Q , (Fig. 3, Panel (b)); this 
is highly consistent with the stellar mass of 10 8 6 M Q found 
by Finkelstein et al. (2010) from best-fitting the observed 
SED. In terms of the halo mass, LAEudf correspond to 
rare, massive 3<r mass fluctuations at z — 8.6 and have halo 
masses ~ 1O 1O,5 M0. The LAEudf have dust masses be- 
tween 1O 5 - 7_6 M as seen from Panel (c) of Fig. 3. These lead 
to values of f c = 0.24 — 0.3 for these galaxies, which translate 
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Figure 3. Summary of the physical properties of LAEs at z = 8.6. In each Panel, points represent galaxies identified as LAEs in our 
simulation; the 7 LAEudf for Xhi = 0.2 (see text for details) are shown by colored symbols. For each LAE, as a function of the observed 
hya luminosity, the Panels represent: (a) the SFR; (b) the stellar mass; (c) the dust mass; (e) the mass-weighted stellar age; (f) the 
mass-weighted stellar metallicity, and (g) the effective ionized region radius around each LAE. 



into Ay ~ 0.75 magnitude, which is in excellent agreement 
with the best fit value of Ay ~ 0.73 magnitude derived 
by Finkelstein et al. (2010) for UDFy-38135539. LAEudf 
have ages between 50-80 Myr (see Panel d), in good agree- 
ment with the value of 10 to a few 100 Myr inferred obser- 
vationally by Lehnert et al. (2010). These ages imply that 
the progenitors of these galaxies started forming as early as 
9.2 < z < 9.7. The stellar metallicity of LAEudf is 0.03- 
0.12 Zq (see Panel e), even at z = 8.6; since heavy elements 
are the main constituents of grains, this is an additional 
argument in support of the dusty nature of these galaxies. 
Indeed dust seems to be required by the measured colors. 

For the 7 LAEudf, ranges between 0.24 - 0.42 
physical Mpc (pMpc), as shown in Panel (f), leading to an 
IGM transmission of T a ~ 20%. If the effects of clustered 
sources is neglected, the nominal individual ionized bubble 
radius is ~ 0.24 pMpc, which makes T a ~ 12%, render- 
ing these galaxies undetectable in the Lycv; the effect of 
source clustering is crucial for making these galaxies vis- 
ible as LAEs. Although the galaxies clustered around the 
7 LAEudf are too faint to be seen in the Lya, nonethe- 
less their UV continuum can be bright enough to be de- 
tectable via standard dropout techniques. In Fig. 4 we show 
the probability (at limiting magnitude Muv ^ —16) of such 
clustered galaxies. Averaged over the 7 LAEudf, there is a 
70% and 15% probability of JWST and HST/WFC3 (with 
limits of Muv <= —16 and <= —18 respectively) finding 
such a clustered galaxy in a radius ~ 0.4 pMpc (or 86 arcsec) 
at z = 8.6. 

As for the contribution of LAEudf to reionization, 
the photon rate density required to balance recombinations, 
g rec , can be expressed as (Madau, Haardt & Rees 1999) 
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Figure 4. Average probability of finding a galaxy at least as 
bright as the magnitude Muv plotted on the x-axis (in bins of 
0.5 magnitude), that shares a common ionized region with any of 
the 7 LAEfjjjF- The error bars show the Poissonian errors. 

9 " ec = 10 5^ ) {0^22) S MPC ' (1) 

where C is the IGM clumping factor. For the 2% HI ioniz- 
ing photon escape fraction used in this work and assuming 
C — 5, the LAEUDF have an HI ionizing photon output 
~ 10 47 ' 2 s _1 Mpc -3 ; at most, such objects could have con- 
tributed ~ 0.01% of the photons needed to balance recom- 
binations at z = 8.6. 

Finally, we estimate the probability of finding a LAE 
like UDFy-38135539 in the experimentally sampled volume 
(1 x 10 4 cMpc 3 ). Our simulated volume (1.08 x 10 6 cMpc 3 ) 
is about 100 times larger; therefore we sample the simulated 
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volume by randomly placing the observed volume within it. 
We find that only 7 of these sub-volumes contain 1 LAE 
similar to UDFy-38135539 for \hi = 0.2. This translates 
into a detection probability of only about 7%, classifying 
the discovery as a relatively serendipitous one. 



4 CONCLUSIONS AND DISCUSSION 

We have constrained the IGM ionization state at z = 8.6. 
We find that without requiring any fine-tuning of the model 
parameters from z = 6.6 and for the maximum possible 
Lya luminosity emerging from the galaxy itself, no LAE are 
found in the observed luminosity range of UDFy-38135539 
for XHI > 0.2. For xhi = 0.2, we find 7 LAEs (designated 
LAEjjdf) whose observed Lya and UV luminosity fall in 
the ranges of UDFy-38135539 (Lehnert et al. 2010) and ID 
125 (Finkelstein et al. 2010). LAEudf are observable in the 
Lya only because an overlapping of their H n regions with 
those of their nearby galaxies make the effective H n radius 
~ 0.24 — 0.42 pMpc; averaged over the 7 LAEudf, there 
is a 70% (15%) probability of such a clustered source being 
found by JWST (HST/WFC3) observations. 

We have also bracketed the physical properties of 
UDFy-38135539: the SFR range between 2.7 - 3.7 M© yr _1 , 
it has 10 8,3-8,7 Mo of stars with a mass weighted age of 
50-80 Myr and a mass weighted stellar metallicity between 
0.03 — 0.12 z?©. It is dust enriched with dust masses ~ 
1O 5,7 ~ 6 M0 which translates into Av ~ 0.75 mag and a color 
excess of E(B — V) ~ 0.16 — 0.19 using the supernova dust 
extinction curve. 

We add three caveats. The first concerns uncertainties 
related to dust attenuation of Lya photons. As has been 
shown by many works (Neufeld 1991; Hansen & Oh 2006; 
Finkelstein et al. 2008; Dayal et al. 2009; Dayal, Ferrara & 
Saro 2010; Dayal, Maselli & Ferrara 2011), the relative es- 
cape fraction of Lya and continuum photons depends on the 
distribution (smooth/clumpy) of dust in the ISM of high-z 
galaxies. Lacking data at z ~ 8 to support either of the 
two possibilities, here we have used the maximum value of 
f a = 1. Secondly, an increase (decrease) in the escape frac- 
tion of H i ionizing photons (/esc) would decrease (increase) 
the intrinsic Lya luminosity, while increasing (decreasing) 
the sizes of the H n regions built by each galaxy, thereby af- 
fecting the observed Lya luminosity. Thirdly, outflows might 
allow Lya photon transmission also from a substantially 
more neutral IGM (Dijkstraet al. 2011) than the upper limit 
of Xhi > 0.2 found here. Although interesting, such result 
has been obtained under the highly idealized conditions of 
a smooth, spherically symmetric expanding shell of HI gas, 
whose radius and column density have been tuned to suit- 
able values. As it is very likely that shell fragmentation oc- 
curs as a result of Raleigh- Taylor and Kelvin-Helmholtz in- 
stabilities arising from the complex gas velocity field around 
the galaxy (e.g. see Fig. 2 of Fangano et al. 2007), such sug- 
gestion needs to await more detailed investigations. 
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